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Abstract 

The result that closed systems evolve toward equilibrium is derived entirely on the ba- 
sis of quantum field theory for a model system, without invoking any of the common extra- 
mathematical notions of particle trajectories, collapse of the wave function, measurement, or 
intrinsically stochastic processes. The equivalent of the Stosszahlansatz of classical statistical 
mechanics is found, and has important differences from the classical version. A novel result of 
the calculation is that interacting many-body systems in the infinite volume limit evolve toward 
diagonal states (states with loss of all phase information) on the the time scale of the interaction 
time. The connection with the onset of off-diagonal phase coherence in Bose condensates is 
discussed. 



1. Introduction 

The controversy over the formulation of the Second Law of thermodynamics in terms of 
statistical mechanics of particles is well known H). Boltzmann presented his H-theorem as a 
way of connecting the deterministic, time-reversible dynamics of particles to the already well 
known, time-irreversible Second Law. This theorem, in Boltzmann's formulation using classical 
mechanics, says that a many-particle system evolves irreversibly toward equilibrium via collision 
processes of the particles which can be calculated using known collision kinetics. Many scientists 
in his day found his arguments unpersuasive, however, partly because they involved notions of 
coarse graining which were ill-defined, and partly because the outcome of irreversibility seemed 
to contradict well-established ideas of time-reversibility such as Poincare cycles. 

Later on, statistical mechanics was given firmer footing on the basis of quantum mechanical 
states, and this approach is used in most modern textbooks [2|. Typically, the notion of entropy 
is introduced in terms of counting states, and the equilibrium distributions of quantum particles 
are deduced on the basis of particle statistics. The basic approach is the same as Boltzmann's, 
using probability and statistics of particles, with the only difference being that the particles have 
well-defined quantum states. This gives the impression that the notion of indivisible particles is 
crucial for understanding statistical mechanics. 



'Email: snoke@pitt.edu, Phone: 412-624-9007, Fax: 412-624-9163 
Preprint submitted to Elsevier 



December 14, 2011 



In this article, we show that by using modern quantum field theory it is possible to deduce 
irreversible behavior for a closed system entirely on the basis of wave mechanics without ever 
invoking the concept of a random ensemble of particles. Starting with a nonequilibrium many- 
body wave function, an equivalent of the Boltzmann scattering integral can be written down, 
by which expectation values of the wave function describing the system can be evolved deter- 
ministically and irreversibly toward equilibrium in the thermodynamic limit. This method has 
already been used to successfully explain various real experiments, but we examine in detail the 
assumptions which are involved in justifying this approach. 

While one may interpret the many-body wave function as giving the probability of observa- 
tion of discrete particles, it is not necessary to do so in order to calculate any of the values we 
calculate here; the equilibration follows entirely from the mathematical structure of the contin- 
uous wave function. This calculation supports the philosophical view that the field of quantum 
field theory is "real." By this we mean that the notion of particles is a subordinate concept of 
the field theory, useful in some circumstances and dispensable in others. In the case studied here 
of equilibration, the particle concept is dispensable. In general, particles are entirely deducible 
from the field theory as eigenstates (resonances) of the field. The many-body wave function is 
therefore a more fundamental physical description of the system. 

A surprising result of the calculation here is that systems generally evolve toward diago- 
nal states (many-body states with no phase coherence between different single-particle quantum 
states) on time scales comparable to the collision time. Thus, the calculation shows a type of 
"collapse" of the wavefunction that comes entirely from the mathematics of quantum field the- 
ory without need for invoking external observers or other, more esoteric concepts. 

Several other works have addressed the evolution of a quantum mechanical system toward 
equilibrium, notably Refs. [3| and |4|. These calculations have typically taken the approach of 
using a "bath" or "environment" which interacts with a smaller subsystem; the bath is assumed 
to already be in equilibrium. In the present work, we treat the entire closed system as one system, 
with the same equations of evolution for all parts, and all of which can be far from equilibrium. 
The essential puzzle we face is that in this case, quantum evolution is unitary and hence the 
von Neumann entropy of a closed system ought to be a constant of the motion. This situation 
has been studied by Gellman and Hartle 0, who used coarse-graining ideas to study the time 
evolution of quantum systems and the cross-over to classical irreversible evolution. Polkovnikov 
|6 1 has introduced the concept of the microscopic "diagonal entropy," which is the von Neumann 
entropy computed ignoring the contributions from the off-diagonal parts of the density matrix. 
He has shown that the diagonal entropy has all the properties expected of the thermodynamic 
entropy. We will present here a derivation of the quantum Boltzmann equation which yields 
irreversible time evolution and show that the quantum coherences represented by simple off- 
diagonal elements of the density matrix do indeed (generically) decay towards zero. As shown in 
Appendix A, an odd result of the Polkovnikov definition of diagonal entropy is that the entropy of 
a boson system near condensation decreases for some period of time. A generalized H-theorem 
can still be derived with a revision of the definition of entropy that also involves only diagonal 
elements, but which has less obvious connection to the von Neumann entropy. 

A key concept in modern quantum information theory that helps us resolve the paradox de- 
scribed above is that of entanglement entropy 00. A strange feature of quantum mechanics 
is that a system can be in a pure quantum state with zero entropy even though each of its sub- 
systems has finite entropy. This is because the entanglement between the different subsystems 
gives non-classical correlations which make a negative contribution to the entropy of the overall 
system. For example, the reduced density matrix for a portion of a quantum system in its ground 
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state obtained by tracing over the other parts of the system describes a mixed state with finite 
energy fluctuations and finite entropy [8|. Coarse graining, in which we trace over short distance 
degrees of freedom to obtain a long-wavelength description [5], can thus lead to the appearance 
of decoherence, classicality, and entropy. 

We discuss below how apparently irreversible decoherence naturally arises within the ap- 
proximations made in deriving the quantum Boltzmann equation. A new result here is that the 
average values of the off-diagonal elements of the generalized density matrix, which hold the in- 
formation about the entanglement, all tend exponentially quickly toward zero. In the low-density 
limit, in which both bosons and fermions recover classical statistics, the ratio of the off-diagonal 
elements to the diagonal elements strictly vanishes. Thus, while the total von Neumann entropy 
does stay constant, the entropy in the off-diagonal terms becomes spread over so many possible 
entanglements that it becomes inaccessible to experiments. 

2. Derivation of the Quantum Boltzmann Equation 

The equation which describes the evolution of a many-body system toward equilibrium 
is commonly known as the quantum Boltzmann equation; to be accurate, the term "quantum 
Boltzmann equation" should be used for an equation which describes both the real-space and 
momentum-space evolution of a many-body system, while the term "Fokker-Planck equation" 
should be used for an equation for the evolution in momentum space only, but in the case of a 
homogeneous gas, these are the same, and the two terms are often used interchangeably. The 
full classical Boltzmann equation is usually assumed to include terms for the spatial drift of the 
particles (see, e.g., Ref. [9|, Section 5.9). Therefore we will use the term "classical Boltzmann 
integral" when referring below to the analogous classical calculation for the evolution of the 
particle distribution function. 

A basic result of quantum field theory which we will use here is that a quantum field can be 
written in terms of second-quantization operators which obey the commutation relation 

[a k ,a\] = a k a[ - a[a k = 1 (1) 

for bosons, and the anticommutation relation 

{a k , a[) = a k a[ + a[a k = 1 (2) 

for fermions. In a homogeneous system, the subscript corresponds to a wave vector k. (We drop 
the vector notation on subscripts everywhere in this paper.) The mathematics of these operators 
goes back to Dirac IfTOll and is presented in many recent textbooks such as Refs. IfTTI and Ifl2l . 
These operators naturally lead to the language of particles, because the operator a! acting on a 
wave state increases the eigenvalue of the operator N k = aia k by 1 and the operator a, decreases 
the eigenvalue of N k by 1, making it natural to treat N k as giving an integer number of particles, 
and a k as the "creation" operator and a k as the "destruction" operator for single particles. More 
fundamentally, the a k and a k operators represent the complex amplitude of a wave of the quantum 
field in mode k, and N k represents its intensity. We will often use the particle language in this 
paper, however, for convenience. Instead of saying "expectation value of the operator N k ' we 
can just as well say "average number of particles in state k;" instead of "low field amplitude" 
we can equivalently say "low particle density," and instead of "nonlinear mode-coupling terms 
in the wave equation" we can refer to "particle-particle interaction terms." 
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2.1. Quantum Boltzmann Evolution for a General State 

A quantum field in the limit of zero nonlinear terms, i.e., zero particle-particle interactions, 
has eigenstates known as Fock states, or "number" states, which are defined by starting with a 
vacuum state and repeatedly using creation operators, as follows: 

i-r (4f 

Here Nk is an integer which gives the occupation number of each state k. (For fermion operators, 
we must apply the restriction Nk < 1.) The "vacuum" state |0) is the zero-particle state, which 
means the ground state of the system, with no excitations. The set of all Fock states makes up a 
complete set of states, which means that any quantum state \ifr) can be written as a superposition 
of Fock states: 

\ifr) = ^ j a n \n), (4) 

n 

where the states \n) are Fock states and a n is the phase factor for each state. 

In a non-interacting system, the total Hamiltonian is given by Hq = 2 Nifiwk + Eo, where 
Eq is the energy of the ground state. If there were no additional terms in the Hamiltonian, 
eigenstates \n) would keep the same number Nk of quanta in each state k for all time, and there 
would be no change of the expectation value of the number operators. Interactions which allow 
the system to approach equilibrium arise from nonlinear terms which multiply creation and de- 
struction operators acting on different fc-states. Such terms arise from anharmonic terms in the 
Hamiltonian — the non-interacting Hamiltonian is equivalent to each £-mode of the system being 
a perfect harmonic oscillator. 

There are all manner of interactions; here, we assume a number-conserving interaction term 
which corresponds to the collision of two particles, i.e., applying destruction operators to two 
states and creation operators to two other states: 

V=\Yj U ^M,kA 4 a l a k2 a h- ( 5 ) 

In this sum, £4 is fixed by momentum conservation, k\ + = £3 + £4, and U^^mm is me 
associated energy of the transition; it is assumed to be symmetric on exchange of k\ with £3 or 
ki with k\. A derivation of this type of term from wave functions in the infinite-volume limit can 
be found in Ref. (9), Section 5.5. 

This four-operator interaction term may seem ponderous, but is actually the simplest in- 
teraction that allows equilibration in a closed system; an elastic scattering term with only two 
operators does not allow the particles to change energy, while a three-operator interaction does 
not conserve the number of particles, and therefore requires keeping track of at least two terms, 
one of which leads to creation of new particles and one of which leads to absorption of existing 
particles. 

The evolution of the system is found using time-dependent perturbation theory (e.g., Ref. IfTZlO . 
If the initial state of the system is and the state of the system at some later time t is \i// t ), the 
change in the average number of particles in state k is given by 



dim = (MNMt) - mfikWi) 
4 



(6) 



= (mw min N k e- iHot!% \m) - mNMi) 



where \ifr(f)) = e iH °' /n \if/ t ) and V(f) = e iHot/nfy e -iH t/n in the stan( jai-d interaction notation. The 
operator N k commutes with Hq, by definition. If it commutes with V, there is no change in the 
occupation numbers over time. 

The exponential of the interaction operator can be written out as the series expansion, 



d(N k ) = {if,i\\l-(l/ih) J^V(t')df + (7) 

x |(l/z7z) dt'[N k , V(t')] + (l/ih) 2 £ dt' dt"[N k , V(t')V(t")] + ...j |^>. 
Let us look first at the lowest-order term of this expansion, which is 



d{N k ) = i f ^<^|[#b<ty')M«) 
«« Jo 



«> 

Jo 

= 4<^*'^llfc>- (8) 

For the particular case of two-body interactions given by Q, we can resolve the commutator 
by using 

[N k ,a k >] = -a k 8 k #, [N k ,a[,] = a\6 k #, (9) 
which are valid for both boson and fermion creation and destruction operators. We find 

[N k ,V]= (10) 
2 ^Zi Ud ± Ue ^ { a k a l, a k2 a h ~ a[a[a h a k ) , 

where £3 — k\ + ki - k, and Ud refers to the direct term and Ue to the exchange term, and the + 
sign is for bosons and the - sign is for fermions. 

The lowest-order term of the expansion is therefore 

d{N k ) = ^Y^( U D ± U E )((fti\aiala k2 a k t\i(i i )-{i(r i \a^^^ (11) 

If \iffj) is a Fock state, this term vanishes exactly; if the two creation operators restore the same 
two states that the destruction operators removed, then the operator has the form N k N kl and is 
equal for both terms, which cancel. More generally, if \\fii) is not a Fock state, (Hi depends on 
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the value of "off-diagonal" elements of the form (if/^a^^a^a^ \\Jri). We will argue in Section j^j 
that these terms are generally negligible. 

We therefore move on to the second-order terms of the expansion Q. The first of these is 



d(N k ) = jpf d 'f dt"mV{t')\N k , WWi) 





, , ... . dfixlfiW^^e-^-t'^ftkMe-^' 1 "^). (12) 
« 2 Jo Jo 

We insert the general definition (|4| of the state |t^,) in terms of Fock states, and insert a sum over 
the projection operators for all Fock states, 2_j \ m )( m \ - 1> to obtain: 

m 

d{N k ) = V JL f dt' f dt" <,a,Xn> /Ho '' /S ye^ Ho( '"'" )/f >«)<m|[^,y]e-'' H »'" /f >) 
= Y \ f dt' f dt" a^a^-^^e^-^HAV\m){M[^kM\n)- 

£fn' n JO Jo 

(13) 

Since V now acts just on Fock states, we can resolve: 

£<n'|V|m><m|[4, V]|«> = ^ J] 2(C7 fl ± U E ) ^N k[ N K (l ± Afr 3 )(l ± Afr 4 ) 

i ^(C/z, ± U E )[y/N kl N k2 (l±N k3 )(l±N k )- y/N k N k3 (l + N kl )(l + N kl ))<«>">■ 



X 2 



(14) 



where |n") is the same as state |n) but with 8 states changed by the 8 creation and destruction 
operators. The numbers N k , N kl , etc., (without hats) are the occupation numbers of the fc-states in 
the Fock state \ri). The sum over states \m) is eliminated by the requirement that |m) be equal to 
\n') with 4 states changed. The factor {n'\n") = 5 n ' A " eliminates the sum over states n' . Here we 
have assumed that all four k's for the transition from \n) to \m) are different; if some k's are the 
same there will be correction terms. However, momentum and energy conservation in the case 
of \n) = \n"), discussed below, will make this impossible except in the case of forward scattering 
(no change of any momenta), in which case such terms make no contribution to change of the 
momentum distribution. 

If the eight changes from \n) to \n") restore all the particle states back to their original states, 
so that \n ) — \n), then this formula is simplified considerably. This forces k\ and k 2 to equal k\ 
and k.2 or vice versa (giving a factor of 2), or k\ and kL to equal £3 and k or vice versa. In the 
same way, ky can equal either k or k\ (giving the direct and exchange terms) if k\ and k' 2 equal 
k\ and k%, or ky can equal either k\ or k.2 if k[ and k' 2 equal £3 and k. This leaves us with 



d(N k ) = V \a„\ 2 Y -L f dt' f dt" 

n 2H Jo Jo 



+E h -E H -E k )fl% e -i(E h +E k2 -E h -E k )t"/h 



x(U D + U E ) 2 \N kl N h (l ± N k3 )(l ± N k ) - N k N h (l ± N k2 )(l ± N kl )\. (15) 



The time-dependent factors can be resolved as 



J> v i(I 



- 1 



dt" e 



sin 2 (wf/2) 



(16) 



We now make two approximations. The first is that the energy states are close enough together 
that there is some range da> over which the energy-dependent factors and the density of states can 
be treated as constant, so that we can treat the sum over k states as an integral over energy. This 
will always be true in the infinite-volume limit, in which case the states of the system are a con- 
tinuum. The second approximation is to assume that the state is changing slowing compared 
to the rate of oscillation of the phase factors. Both of these approximations are the standard ones 
for deriving Fermi's golden rule for transition rates, and can be called the Markovian limit; we 
will discuss the breakdown of this approximation below. In this case, the time-dependent factor 
can be written as: 



lim 

t— >oo 



sin 2 (wf/2) 



of 



5{<jj)2nt. 



Taking t as a small time interval dt, we then have for ( 13 1: 



dt 



{N k ) 



2n 1 

~h2 



Y/Vd ± U E ) 2 6(E kl + E k2 - E h - E k ) 



(17) 



(18) 



kite 



A{<p i \N h N k2 {\±N h ){\±N k )\ilf i ) - (iff i \N k N h (l±N h )(l±N kl )\i/f i )) 



This is nearly the standard quantum Boltzmann equation, or Fokker-Planck equation, for two- 
body interactions. The only remaining step is to factor the expectation values to have products 
like (N k )(N k >) instead of (N k N k >) etc. We will address this factoring in the next section. 



In going from ( 14 1 to ( 15 1, we restricted ourselves to terms in the sums where the second 



interaction term completely undid the changes in the state \n') made by the first interaction term. 
If we had not done that, for a given set of k's we would have terms of the form 

2 «>» ( ^N k ,N k > 2 (l ±Afr 3 )(l +N K ) y/N kl N k2 (l±N k3 Xl±N k ) ) l 

n 

a\, a\, a k >a K a\a\a k2 a k , (a„\n)) 



(il/i\a\, a\, avaval ala k ,a k \^i) 



*k' "k. 



(19) 



where the eight operators are those which take \n) to \n"). This is another, fourth-order "off- 
diagonal" phase term like the four-operator, second-order off-diagonal term found in ( fTT] i. As 
with that term, we will argue in Section[5]that these fourth-order off-diagonal terms are negligi- 
ble. 

Last, we must examine the other second-order term in Q. This is 



~ X dt ' L dt " { ^ k ' vvwv 



Mi) 
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= ~i X dt ' L dt " § (<"' l ^ i>(f ' )<>(f " )|n> - (n'\V(t')V(f')N k \n)) 
~^L dt ' L ^ ? a "" a " (N '«' ~ N k){n"\V(t')V(t")\n). (20) 

If we consider terms in which the state \n") equals \n), then this term vanishes, since N',', which is 
the value of Nk in the state \n"), is then equal to Nk- If we consider other terms, then we must deal 
with the same fourth-order off diagonal terms discussed above, which we will argue in Section[5] 
are negligible. 

2.2. Factoring the Number Averages 



Equation ( 18 1 gives the rate of change of any average occupation number {Nk) in terms of 
instantaneous averages of the form {^i^N^Nk^d and {^N^Ni^Nk^i/fi). (The two terms with 
four number operators actually cancel out.) What we would like, however, is an equation in 
terms of only {ift^N^ \iffi), {ifti\Nk 2 \<fri), etc., so that we only need to keep track of the average 
occupation numbers, not any correlations between the numbers in different momentum states. 
This requires us to justify using {t//i\N kl N k2 \i/fi) = (<ffi\N kl \t/f i )(ifr i \N k2 \i//i). 

Since the interactions lead to real-space correlations or anticorrelations, this factorization 
will not, in general, be exactly true, but can be taken as a good approximation in an incoherent 
system. One way to justify this is to look at the implications of reasonable assumptions about 
the spatial correlations. We write the density-density correlation function, 

G (2 \r) = <pWp(0)>, (21) 

where 

r>(r) = i/AFWrt = 

kite 

is the local density operator in terms of the spatial field operators. Substituting in this definition, 
we have 



(23) 



V 2 



As we did in the previous section, we will assume that "off-diagonal" terms {a^a^a^a^), with 

all four k's different, are nearly zero, so that only terms with two sets of k's the same, to make 
number operators, will contribute to this sum. This gives us 

G{2)(f) = 4 Z ( ±e_i * ^ + 1 ) {f/k < (24) 

kite 

where the + sign is for bosons and the - sign is for fermions. 

If there is no correlation in the many-body wave function at separation r, we expect that 
G (2) (r) will have the value 

{ifrUr)ifr(r)){f(0)ifr(0)) = ^ £ ^"'HV ZH^> 



V 2 

p\ (25) 



where p - N/V is the average number density, and we again assume vanishing off-diagonal 
terms which make only the number operators relevant. Integrating over all space, we obtain 

fd 3 r(G ( - 2 \d-p 2 ) = ± I^2) + I^(<^ i ^ 2 >-<^ i X^ 2 >). 

ki kite 

The first sum on the right side can be taken as negligible compared to the second, since it consists 
of just the ki = k\ terms of the first term in the second sum. 

Suppose that G (2 \f) deviates from p 2 only for some limited range in space, e.g., 

G (2) (7)=p 2 (l- e - r2/fl2 ), (26) 

where a is some characteristic length. Then the integral over all space is 

'd 3 r(G^-p 2 )~-p 2 a 3 . (27) 



Then we have 



V 

kite 



(28) 



or 



2(<#*A>-<#U<#fc>) ~ ~Yj(N kl )(N h ). (29) 



V 

k,te hte 



The sum on the right-hand side is small compared to the second term in the sum on the left-hand 
side, as long as a is not extensive with the whole volume. We therefore have 

Yj (<#*. > - (Nk, )(N h )) * 0. (30) 
kite 

In other words, on average, the deviation of {Nkfiki) from (Nk l )(Nk 2 ) is zero. 

We can make the same type of argument for the factorization of (N^N^N^) using higher- 



order correlation functions. (The terms with four iv* operators in (18i cancel.) In general, a 
localized correlation or anticorrelation in real space will lead to no significant £-space num- 
ber correlations. The case of a Bose-Einstein condensate, which has infinite-range correlations, 
requires re-examination of the factoring of the number averages. We discuss some aspects of 
Bose-Einstein condensation below, in Section [B~3l 
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3. Solution of the Quantum Boltzmann Equation 



We now have a tractable equation for numerical simulations. We can simply define the occu- 
pation number (N k ) for all k and calculate the rate of change for each (N k ) at every point in time. 
In this section we examine the numerical solution method in detail. 

In all of the above, we have assumed that the fc-states are defined for a set of discrete k 
values. In the limit of infinite volume, however, the spacing between states k disappears, and the 
distribution function N k is a continuous function defined for all k. In that case, we can set (N k ) 
equal to its average value in a small element of phase space d 3 k. This allows us to convert the 
double summation to a double integral, 



d(N k ) _ 2n 
dt ~ h 1(2. 



—V- f 

2ny] 2 J 



d 3 h d 3 k 2 \U D ± U E \ 2 5{E kl + E k2 -E k - £V) 



X [(N kl ){N k2 )(l ± (N k ))(l ± (N k ,)) - (N k )(N k ,)(l ± (N kl W ± <#fe»] . (31) 

where V is the volume, and the + signs are for bosons and the - signs for fermions. If we further 
assume that the system is isotropic, i.e. the distribution function depends only on the magnitude 



of k and not on the direction, then the angles in the integral ( 3 1 1 can be eliminated by analytical 
integration, to reduce this integral to just a double integral over two energies, which can then be 
computed numerically. We can then write the distribution function simply as N{E k ) = (N k ). At 
a given point in time, the change dN(E)/dt can be computed for all E, and then each N(E) is 
updated according to 

dN(E) 

N(E) -> N(E) + — )—dt, (32) 
dt 

where dt is some small time interval, chosen such that the change of N(E) is small during any 
one update. 

In other words, we have a prescription for determining the approach of a nonequilibrium sys- 
tem toward equilibrium deterministically . Figure [T] gives the solution generated by a computer 
for bosons at low density starting from a nonequilibrium distribution. After five times the charac- 
teristic scattering time r, the distribution does not change significantly — the system has come to 
equilibrium. This is not a Monte Carlo simulation. In this calculation, no random numbers were 
used, and no individual particles were tracked. The distribution function N(E) was evolved de- 
terministically using the isotropic quantum Boltzmann equation, such that each successive state 
was completely determined by the previous state. 

We can easily determine the equilibrium distribution by finding the steady-state solutions of 
the quantum Boltzmann equation. In equilibrium, we must have dN{E k )/dt = for all k. This 
will occur if for all the different scattering processes, the forward and backward rates for any 
process are the same: 

N(E k3 )N(Ek 2 )(l ± N(E kl ))(1 + N(E k )) (33) 
-N(E k )N(E kl )(l±N(E kl )(l ±N(E k} ) = 0. 

This is the principle of detailed balance. It can easily be verified that a distribution of the form 

W = ^WTT (34) 
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Figure 1 : Nonequilibrium evolution of an isotropic gas of bosons at low density interacting with the potential J5J, found 
by numerically solving Eq. (31) . The initial state (a) corresponds to equal occupation of all states up to a cutoff (N(E) is 
proportional to the density of states !D(E) oc yE). The distribution (b) is after each particle has undergone one scattering 
event, on average, i.e., the total amount of change of the distribution N(E) is equal to the total integral of N(E). The 
distribution (c) is after two scattering events per particle. The distribution in (d), after four scattering events per particle, 
is less than 1% different from a Maxwell-Boltzmann distribution. The horizontal energy axis is given in units of kgTa, 
where Tq is the final temperature. From Ref. 1 13]. 



satisfies this condition, where the - sign is for bosons and the + sign for fermions. The constants 
a and /3 are determined by the conditions 



N 



Y J N k = ( N(E)D(E)dE, U = ^j E k N k = J EN(E)D(E)dE, (35) 



where O(E) is the density of ^-states at energy E. These conditions follow from number and 
energy conservation, which follow from the fact that the total number operator jV and the total 
energy operator U commute with the Hamiltonian. We can find a and /? in terms of standard 
thermodynamic quantities by noting that if a » 1, the equilibrium distribution is equal to N(E) — 
e~ a e~P E . Equating this to the standard equilibrium distribution N(E) = e 1 *^" 7 e~ E l k " T , which is 
called the Maxwell-Boltzmann distribution, we then have /3 = 1 /ksT and a = -fi/lcsT, or 

W - e C*-JLr Tl . W 
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which is the well-known equilibrium distribution for quantum particles; when the sign is negative 
it is called the Bose-Einstein distribution and when the sign is positive it is called the Fermi-Dirac 
distribution. The quantum Boltzmann equation is therefore another way of deducing the standard 
equilibrium distributions. 

The Maxwell-Boltzmann limit fj. <zz Eq corresponds to low particle density both the Fermi- 
Dirac and the Bose-Einstein distributions. The equilibrium solution shown in Figure [T] for the 
distribution of the particles at late times converges to the Maxwell-Boltzmann distribution times 
the density of states of the particles, which for the three-dimensional gas in the model is propor- 
tional to y[E (see, e.g., Ref. Section 1.8). 

It is easy to show that the equilibrium solution of the quantum Boltzmann equation is the 
same even if we invoke some process other than the two-body elastic scattering mechanism 
assumed here. For example, suppose we have an electron-photon interaction of the form 

= 2l i {al+a p )al 2 a e i, (37) 

where a e \ and a e2 are fermion destruction operators for electrons in states 1 and 2, respectively, 
and a p is the operator for a photon which couples the two states. This gives the detailed balance 
equation 

iV e i(l - N e2 )N p -N e2 {\-N e {){\+N p ) = 0. (38) 

If the electrons are in a Fermi-Dirac equilibrium distribution, then the ratio of the electron popu- 
lation terms is 

N e2 (l-N e i) = e - iE2 -E0/k B T (39) 

N el (l-N e2 ) 

which then implies 

which is the Planck distribution, equal to the Bose-Einstein equilibrium distribution with a = 
0. Therefore, just as with any other equilibrium distribution, we can also derive the Planck 
distribution entirely from the field theory. Although the Planck distribution historically played a 
major role in convincing people of the particle nature of waves, we see here that it fundamentally 
arises because of the field operator commutation relation ([TJ, which led to the final states factor 
(1 + Nk) in ( 18 i. The commutation relation, in turn, follows from the wave nature of the quantum 
field, without reference to any particle collapse or measurement (see Ref. [9 1, Sections 4. 1-4.3). 

The nonequilibrium evolution predicted by the quantum Boltzmann equation has been veri- 
fied experimentally 1151 . Figure [2] shows an example of the evolution of a population of excitons 
in a semiconductor. Modern ultrafast optics experiments can resolve the energy distribution of 
these particles on time scales short compared to the time to equilibrate. Therefore we can observe 
the approach to equilibrium and compare it to the evolution given by the quantum Boltzmann 
equation. As seen in this figure, there is a very good fit of the experiment to the theory. Only one 
fit parameter was used to fit the data at all times, namely, the matrix element for the strength of the 
interaction. The initial distribution was determined by the laser pulse which created the nonequi- 
librium distribution. Long after the laser pulse, the particles remain in a Maxwell-Boltzmann 
distribution. 



4. Where Did Irreversibility Enter In? 

The universal evolution toward equilibrium seen in the numerical solutions of the quantum 
Boltzmann equation is equivalent to the content of the H-theorem of Boltzmann, which says 

12 



Photon energy (eV) 



Figure 2: Solid lines: energy distribution of excitons in the semiconductor CU2O measured at various times following 
a laser pulse with temporal width 2 ps and maximum intensity at time t = 0. Dashed lines: solution of the quantum 
Boltzmann equation for the time evolution of the population using a deformation-potential matrix elemennt for exciton- 
phonon scattering. The theory gives a Maxwell-Boltzmann distribution at all late times. (From Ref. ( 14 1). 
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that when entropy is calculated for a nonequilibrium system, it always increases until the system 
reaches equilibrium. Appendix A shows the connection of the quantum Boltzmann equation 
to the H-theorem in terms of various definitions of entropy. As seen in the previous section, 
though, the basic result of irreversible equilibration does not depend on any particular definition 
of entropy. 

But in a closed, energy-conserving system, we expect the dynamics will be time-reversible. 
In other words, if start with state |/) and evolve to a state \j) at some later time, then we expect 
that if we started with the system in state \ and reverse all the momenta, we should get back to 
state 

The results of the previous sections appear to contradict this view. We began with a closed, 
energy-conserving system, namely a quantum mechanical field Hamiltonian with no interactions 
with any external system, and deduced the expectation values of the many-body wave function as 
it evolves deterministically toward equilibrium according to the proper wave equation; we never 
invoked collapse, measurement, observation, or randomness. From beginning to end we have 
treated only the wave function, without invoking particles at all except to identify them as the 
natural energy eigenstates of the system. Yet we get an irreversible, deterministic approach to 
equilibrium. How have we gone from a time-reversible system to one with apparent irreversibil- 
ity? 

Of course, along the way we have made several approximations in order to generate the 
quantum Boltzmann equation. Is it possible that these approximations are responsible for the 
irreversibility? Let us now examine in detail the assumptions which were made. 

Let us first address an apparent inconsistency which may already have occured to the reader. 



On one hand, in using ( 17 1, we assumed that the time interval for the evolution was long enough 



that the t — > oo limit could be used, while in writing ( 18 1 we assumed that t is short. Both of these 
limits can be valid. We pick dt short enough that N(Et) of any given state is almost unchanged, 
but also long compared to a cycle of the phase factors which are relevant in that energy range. 

This natural energy range is determined by the variation of the occupation number N(Ek) 
with energy. If N{E k ) is a continuous function of k, then there will always be some energy 
range (Ek,Ek + dE) over which N{Ek) can be considered constant. This defines the "natu- 
ral" energy range AE, which in turn defines a natural frequency to = ti/AE. For the function 



sin (a>t/2)/(2naj t) to be well approximated by the delta-function 6(a>) as we did in (17), we 
therefore want to choose dt such that 

1 dN(E k )Y l , H 



, <K dt <K — . (41) 
N(E k ) dt ) AE 

It is possible to have physical situations, however, in which the rate dN(Ek)ldt is large, 
forcing the time interval dt to be small, while N(Ek) is strongly varying over a narrow range of 
energy, so that the natural energy range AE is small, and therefore the period 1 /a> is large, such 
that there is no range of time which satisfies both assumptions above. In this case one cannot use 
the quantum Boltzmann equation. Instead, one must use a non-Markovian evolution equation 
with a quantum memory kernel, in which memory is included of the state of the system at earlier 
times. Such methods, generally known as "quantum kinetics," have been developed at length by 
Haug and Koch [16|. The general effect of using quantum kinetics is to introduce oscillations in 
the response of the system on very short time scales, typically of the order of femtoseconds for 
solid state systems. The solutions always approach an equilibrium distribution irreversibly on 
long time scales, however. This is because there will always be some scattering into regions of 
phase space where N(Ek) is slowly varying enough for the quantum Boltzmann limit to apply. 



Another approximation which we used in deriving the quantum Boltzmann equation was 
restricting the perturbation theory to the leading order in the series expansion (|7J. Would in- 
cluding higher-order terms remove the irreversibility? No, these higher-order terms just give 
multi-particle interactions instead of two-body interactions. Higher-order expansions of Fermi's 
golden rule are well known (see, e.g. Ref. [9 |, Section 8.1) and just give higher-order corrections 
to the matrix element used in the transition rates, i.e. corrections to the scattering cross section 
beyond the Born approximation. 

Some may object that the quantum Boltzmann equation only keeps track of the average value 
of N k and does not keep track of fluctuations. There are two reasons why accounting for fluctu- 



ations will not restore reversibility. First, as discussed in Section 2.2 we expect the numbers of 
particles in different fc-states to be uncorrelated, which means that in an element of phase space 
(fk which includes many fc-states, the number fluctuations should be governed by Poisson statis- 
tics for uncorrelated systems. In this case the number fluctuations are proportional to -\fN\Etj, 
and therefore the ratio of fluctuations to the total number decreases as 1 / -\fN\Ekj, making them 
negligible in any macroscopic system. Furthermore, fluctuations do not lead to periodic behav- 
ior. The mean distribution N(E k ) will still obey the evolution given by the quantum Boltzmann 
equation. 

Where did the irreversibility enter in, then? The place to look is the iteration procedure 



(32 1. Note that the quantum Boltzmann equation ( 18 1 gives the change in the average occupation 



number Nk, using only the information of the average occupation numbers for the states. In the 
derivation of the quantum Boltzmann equation, we threw out terms of form 



(2) 



= {\l/i\a\ala ky a k \^i), 



and 



One way of speaking of this is to say that we keep only the "diagonal" information of the state of 
the system, and discard the "off-diagonal" information. This terminology is a generalization of 
the density matrix formalism, where the first-order density matrix of a quantum system is defined 
as (see Ref. |9 |, Section 9.2) 



(42) 



(43) 



P = 



( {a h a k) 
(a\a k ) (ala k2 ) <a*a fe > 

(a\a,) (4a,) (a'U k ) 



(44) 



The diagonal elements of this matrix are the average occupation numbers, while the off-diagonal 
terms account for phase correlations in the system. We can generalize to define as "off-diagonal" 
any expectation values which do not involve only matched pairs of creation and destruction 
operators of the form a} k a k . 

Suppose that the system starts in a purely diagonal state, with all off-diagonal terms of all 
orders strictly equal to zero. After a short time dt, while the quantum Boltzmann equation gives 
the correct value for the change of (Nk), the state of the system will have evolved to a different 
superposition of Fock states, and therefore the off-diagonal terms may be nonzero. The procedure 
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( 32 1 amounts to replacing the state of the system with a purely diagonal state after every time 
step dt. 

Let us compare this to the standard derivation of the H-theorem in terms of classical particles. 



The classical Boltzmann integral looks just like the quantum Boltzmann equation (31 1 derived 



above, in the limit when the particle density is low so that the (1 ± Nk) final states terms can 
be dropped. Just as we have used a diagonal state with a continuous (N^) for the momentum 
distribution of the particles, the classical Boltzmann integral describes the state of the system 
in terms of an continuous velocity distribution function, which keeps track of the number of 
particles with a given velocity, but not the exact position and velocity of individual particles. 

Suppose that at time t — we have an isotropic, nonequilibrium velocity distribution of 
classical particles, and we evolve this distribution to a new velocity distribution at a later time 
time t = dt via the classical Boltzmann integral. Following the same iterative approach described 
above for the evolution of a system using the quantum Boltzmann equation, we can then use the 
classical Boltzmann integral to evolve the distribution to the next time step, and so on until the 
system reaches equilibrium. This was Boltzmann's approach. But since the microscopic laws 
are time-reversible, a problem arises. If the velocities of all the particles at time dt are reversed, 
then the system must evolve back to its original state at time t — and away from equilibrium. 
The classical Boltzmann integral will never show this behavior, however, because the overall 
velocity distribution for an isotropic system would be unchanged by reversing all the individual 
velocities. 

The resolution of this which gained acceptance in statistical mechanics is the Stosszahlansatz 
(for a review of this history, see Ref. HI). This assumption says that while some states of the 
system do exist which would give time-reversed behavior, for any given velocity distribution the 
number of such states is extremely small compared to the total number of states with the same 
average velocity distribution. Effectively, then, at each time step in the iteration of the classical 
Boltzmann integral, we replace the exact velocity distribution of the particles with an average 
distribution. This loss of information ensures the irreversibility of the evolution ifTTl . In this 
approach, the H-theorem is thus seen as a statement of high probability, not absolute necessity. 
The chances of having a set of velocities which exactly reverse the time dynamics and move the 
system away from equilibrium is vanishingly small, but not strictly zero. 

The same loss of information happens in the numerical evolution of the quantum Boltzmann 
equation. Keeping track of all the phase factors of the superpositions of different Fock states in 
the full many-body quantum state is like keeping track of all the individual particle velocities 
in the classical case. If we did keep track of these terms, then we could write down the exact 
quantum mechanical state at time dt, and using this as an initial state, by reversing time we 
could recover the initial state at time t — 0. Replacing the state of the system with a diagonal 
state at each time step in the evolution of the quantum Boltzmann equation is like replacing the 
exact velocity distribution with an average velocity distribution at each time step when using the 
classical Boltzmann equation. The loss of the off-diagonal information, known as "dephasing," 
ensures irreversibility. 



5. Calculation of the Off-Diagonal Time Evolution 

As discussed above, the elimination of off-diagonal phase information plays the same role 
in the quantum calculation as the Stosszahlansatz in the classical calculation. The quantum 
calculation allows to us go one step further, however. Since we have a closed system with the 
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full quantum mechanical many-body wave function written down, we can compute the dephasing 
rate. Then we can ask how valid an approximation it is to discard the phase information as we 



did in Section 2.1 In that section, the evolution of (Nk) depends, in principle, on the values of 
all of the different p (2) and p (4) terms. We note that the evolution of (Nt) depends only on these 
terms, that is, on averages of the phase factors, and not on detailed knowledge of single values 
of the phase factors. 

In the rest of this section, we show how to do this. We will focus on the evolution of the 
p (2) terms, and return to discuss the p (4) and other terms at the end of this section. We define the 
operator pf}, .„ ,,„ = a k a k ,a k „a k ,„, and assume that the four momenta are all different, because the 
case when two or more momenta are the same will be greatly suppressed due to the requirement 
of momentum conservation in the interaction Hamiltonian which gives rise to evolution of the 
p (2) terms. 



Following a similar approach to that of Section 2. 1 we write for a given off-diagonal term, 



mt)\pfl r , k -At)\m) - mpfl rr ^i>- (45) 



As we did for (|7]), we expand the time-dependent exponential operators as 

<^fi'i"i-> = m (1 - l - h f o nt'w + -L JT* dt> JT dt"v(t')nt") + . . )p% rr xt) 

' 1 + - f V(t')dt' + f dt' f dt"V(t')V(t") + 
\ in Jo (in) 1 J,, Jo 



jo \ in r jo jo 1 

-(^fl^M^i)- ( 46 ) 



The lowest-order term is 



the first-order terms are 



<M(PTu M ,„M-pZ», M ,,W i )\ (47) 



J k,k',k",k 

JO 

while the second-order terms are 



4 f dt '^>i 

in Jo 



\[pfl k „ k ,M)Mt')M,)\ (48) 



\ J df I dt " & A ^ f &wr#»®V(f'Mi) 



n 2 

1 



dt' dt"(fr\W)W)pfi k „ Ami) 

Jo 



(49) 



We cannot assume that jx. k „ k ,„ commutes with Ho, as we did for Nk in the derivation of the 
quantum Boltzmann equation, so we must keep all of the e lH "'^ n terms in the definitions of the 
time-dependent operators for the moment. In the following sections we compute each order of 



the expansion (46 1 separately. 
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5.1. Zero-Order Off-Diagonal Evolution 

As before, we write \tj/i) as a sum of Fock states, to have the zero-order term 

n,n' 

(50) 

Applying the time-varying exponential factors to the Fock states gives us 

- + £jf - £*" - £*»') £ a;,a„<n'|p^, >jt „ iit ,„|n> 

= + £*, - E r - ^"X^^^l^}. (51) 

Thus, if there are nonzero p (2) terms, these will rotate in phase proportional to the degree of 
violation of energy conservation. 

5.2. First-Order Off-Diagonal Evolution 

We start by working on the first-order term. Inserting the definition of V, we have 

wfl^xt),W)m) = (52) 

k u k 2 ,k 3 

where we abbreviate the momentum-dependent matrix element simply as U. As we did in deriv- 
ing the quantum Boltzmann equation, we write the initial state as a superposition of Fock states 
of the form Q: 

£ a* n ,a n J] |(<n> ! ' ffo ' /fi aI4a^^^ e - iff0,/s e iff ^ /s 4 i 4 2 a t3 a it4e -^'' s |n> 

n,n' ki,k 2 ,kj 

-<n>^^ fi 4 ia t ja% ^ e -iHo^s e Wo ? /» a t a t a ^,^„ e -iHo«/S| n>) . (53) 
For n' = n, this becomes 

£ \a/2(U D ± U E ) (^E v -Ey,-E v ,,) it -ni%^n)^n) {l ± ± ^n)) 

_ e -/(£ i+ £ t ,-£ t ,,-£ 1 ,,0(/-O/^W A? («) (1 ± ± jy-M)) ( (54) 

where, as in Section [2~T] Up and Ue are the direct and exchange interaction constants, and the 
numbers Nj* , , etc., without hats, are the occupation numbers of the ^-states in the Fock 
state \n). The ± signs here give the cases of bosons and fermions, as in the quantum Boltzmann 
equation. If n' + n, we are left with off-diagonal terms, either p (T> or higher-order phase terms, 
namely a six-operator p (3) term or an eight-operator p (4) phase term. We will discuss these below. 
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Performing a time integral for the case n' = n, we have 



fdt'e 
Jo 



1 - I0)t + . . . - 1 
(1 + /wf + . . .) =a t. (55) 



Thus, for a small time interval f, the first-order evolution from a diagonal state gives 



d_ 

dt 

2i 



(U D ± U E )m (N k „NMl ± N k )(l ± N k >) - N k N k ,(l ± N k „)(l ± #*»)|fc>) • 



This term gives us phase coherence that accumulates proportional to the net scattering rate con- 
necting the four states, but without any condition of energy conservation. This term is imaginary, 



which multiplies the imaginary prefactor in ( 1 1 1, giving a real change of the expectation value 



(N k ). Since this term increases linearly with t, if the phase factors start out at zero value, the 
change in population in ( 1 1 1 will increase as f 2 , which will be negligible on short time scales. 



But if the phase factors grow in time to have nonzero values, then they could contribute to the 
evolution of (Nk). 



As mentioned above, if we do not require \n') = \n) in equation (53 I, we will have terms with 
higher-order off-diagonal terms instead of just number operators N k . In this case we will have 
terms like the following: 



(n>^' fl 14«,„ % „ e -^ f e ^''«l 4 4 3% «, ie -^ f, |n) 
-(n'\e iHot ' a\a\ a h a k e~ im ' e iH °' 'alal,a k „a km e~ mt \n). 

(57) 



If none of the operators a kj , , a£ , or act on the same ^-states as the momenta in p^, k „ k ,„ , 

i.e., if all the operators in the interaction term commute with pfj. k „ k ,„, then these two terms will 
cancel. 

(2) 

To deal with the case when the interaction term does not commute with p, L ,„ ,,„, we must 
deal with a number of possibilities. Appendix B treats these in detail. The general result is that 
if there are nonzero p {1) factors, then these can contribute to the first-order evolution of p[ 2 ], ,„ ,„, . 
These other p <2) off-diagonal factors always appear in sums, such as 

J]f(N k ,N k ,,N e ,,N^)pfl MM , (58) 

h,k 2 

where f(N k ,N k ',N k »,N k '») is some function of the occupation numbers, which can be either 
positive or negative. We might in general expect that these sums will average to zero, but as we 
will see below, even if they do not, there will be an overall evolution term to suppress these p (2) 
factors. 

5.3. Second-Order Off-Diagonal Evolution 

We now proceed to the second-order terms (49 1. First, we examine the time integrals. We 



assume that the operators in the first V(t') undo whatever changes have been done by V(t"); 

(2) 

otherwise we will not be left with Wvp k Lu,, k ,„ wu- Terms which do not do this will involve 
higher-order off-diagonal phase correlation terms pP\ p (4) , etc. These terms in principle can 
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contribute to a second-order contribution to the evolution of (i^, |p (2) |i/'/), but as with the first-order 
calculation, we take a perturbative approach that the p (2) terms will dominate over higher-order 
off-diagonal terms. What we are interested in now is any terms proportional to (t^,|p (2) |t/';), so 
that we can see what happens to existing off-diagonal phase accumulation. 

If V(t') undoes whatever changes have been done to by V(t"), then each of the three 
terms in ^ gives a factor e KE h +E h -E h -E h )t' /n Xe -i(E h ±E h -E h -E ti W '/ft j but in the first case f and 

t" are integrated from to t separately, while in the last two terms t' is integrated from to t but 
t" is integrated from to f . The first case gives a factor (2nt/fi)6(E kl + - E k , - E ki ) as given 
in ( fT7| ). The second case can be written as 

f dt' f dt" e ,w V w " = f dt' f dt" e'"''e~ ia "" - f dt' f dt" e^e' 10 "" 
Jo Jo Jo Jo Jo J/' 

Jr*t r*t r*t r*t—f 

dt' \ dt" e UMf - , " ) - I dt' I dt'" e 
o Jo Jo Jo 

= f dt' f dt" e ia "'e- ia "" - f dt"" f 
Jo Jo Jo Jo 



-Uo(t—f)i6)f 



dt'" e- m e lh " . 

(59) 

The second term on the right in the last line is equal to the term on the left, and therefore we have 



f dt' f dt" e ia "'e- ia "" = i f rff' f rff"e ,w < 
Jo Jo 2 Jo Jo 



7rf 

-r^iEti + E k - E k} - E ki ), 
n 



(60) 



i.e., 1 /2 of the first case. 

For convenience we write one term of the V sum as V, = Uat a k a^a^ . After removing the 
time dependence, the second order term of interest is 

m (%P i2) Vi - \%Vfi 2) - \P (2) %V) l<A;>. (61) 

This is the generalized Lindblad operator for a many-body system, used often in quantum optics 
|[T8l . If Vi and % both commute with p (2) , then this vanishes. Thus, all of the terms in V which do 
not include a creation or destruction operator with at least one of the four momenta k, k',k", k'" 
in p (2) do not contribute to any change of this dephasing term. This makes sense, since only 
interactions with these states should contribute to this dephasing term. 

~ (2) 

We must therefore consider each type of V, term that does not commute with p k j, k „ ,,„ . We 

(2) 

will take these in order: terms with one momentum the same as in P k k , k „ k ,„, terms with two 

(2) 

momenta the same as in p k k , ,„ k ,„, terms with three the same, and terms with four the same. 

~ (2) ~ 

The interaction term V, will not commute with p k L k „ k „, if one of the four operators in V, 
acts on one of the same momenta as the four in p ( kk , k „ ,„, . We assume that has the same four 
operators as in V,, so that all the changes in |^,-) made by V, are reversed, as discussed above. In 



this case, (61 1 will have terms like the following: 

1 + t ! , ! 
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; - ^44^--^, (62) 



which is true for both bosons and fermions, using the fermion property that a double creation 



operator vanishes. The remaining operator in (62 1 goes back into the definition of p (2) , leaving 
behind three operators in each of V, and V,-< . 

Setting k\ — k in %, we can pick the momenta in % four different ways to restore the same 
Fock state: 

^A^^kAA^^k + a A kl a h a kAkAh a k 2 a k 

+a\a\a k a k a\a\a k a k + a\a\a k a k a\a\a k a k . 

(63) 

If U is a constant, this just gives a factor of 4 multiplying the (U/2) 2 prefactor that comes from 
the two V terms in second order; more generally it gives the prefactor | Ud(Ud ± Ue)- There will 
be another four terms giving ±\Ue(Ud ± Ue) if we take ki — k, thus giving (Ud ± Ue) 2 , which 



is the same matrix element as in the quantum Boltzmann equation ( 18 1 



Setting k-$ or k\ in V,- equal to k gives two more sets of four terms, using 

1 t t 1 , + 1 



a k a\a\ - -a k a\a\ - -a\a k a\ = ±-a\. (64) 



where the + is for bosons and the - is for fermions. 

(2) 

Thus, for the momentum k in p\ L ,„ k ,„, we have a factor of the form 

l -(U D ± U E ) 2 [±#*A(1 ± N h ) - N k2 (l ± N k3 )(l ± N h )] . 

Since k\,kz,kj,, and k<\ are dummy variables, we can switch k\ with k.2 and k^, with k\, without 
changing the matrix element. We thus finally have, for the second-order term, 



J t (P ( kkk",k^ = Y (Ud * UeY 2 2j 6(Ek + Ek * ~ Ek > ~ Ek < } (65) 

x <^\pfkk",k'" [±^3^ 4 (i ± K) - Ko. ± ^xi ± no] iiA,-> + . . . . 

where the . . . indicates that there are three more terms of exactly the same form, for k', k", and 



In 



(U D + U E ) 2 \ V 6(E k + E h - E h - E h ) 



;/// • -(2) 

k in " K ' 



Pk,k',k",k" 



Under the same assumption of no correlations that we used to factor ( 18 I, discussed in Sec- 
tion ! 



2.2 we can factorize (65 1 into products of expectation values. We then have 



J t <Pflr,k-> - { P Z rri ) 2 -^(U D ±U E ) 2 \Y J ^E k+ E kl -E h -E k4 ) (66) 

ki,h 

x [±<^,)<^ 4 >(1 ± (N k2 )) - {N k2 )(l ± <iV fc3 »(l ± <^ 4 »] + . . . . 

The population factor {N k} ){N ki )(l +{N k2 ))-{N k2 )(l +{N k3 ))(l +{N ki )) does not vanish in equilib- 
rium, unlike the comparable term in the quantum Boltzmann equation, because it does not have 
the balancing factors (N k ) and (1 ± (N k )). For fermions, this rate is always negative, leading to 
dephasing, while for bosons, it can be positive if the net scattering rate into one of the ^-states is 
positive. Net scattering into a state can occur, for example, when there is macroscopic coherence 
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due to Bose condensation, which gives a macroscopic number of particles in a single state. Ap- 
pendix C discusses the case of Bose condensation. At low density or high temperature, however, 



when (Nk) <s 1, the out-scattering term in (65 1 will dominate over the in-scattering term, since 



the in-scattering involves a product of two occupation numbers, while the out-scattering term has 
only a single occupation number. In other words, at low occupation per state, given a particle in 
state k, it is much more likely for that particle to scatter out of that state than for another to enter, 
in a given time interval. Therefore at low density the overall rate is always negative. 

5.4. Summary of Dephasing Calculation Results 



We have seen that on the one hand, that (56 1 implies that scattering does lead to deviation 



from a pure diagonal state, with an increase of the phase factor (p (2) ) proportional to the net 
scattering rate. But in the case of fermions, or bosons at low density, that phase coherence can 
never build up; the result ( |66] > is of the form 

>>> = -^>, (67) 
at t 

which gives exponential decay of the phase with a time constant r; the value of r is equal to the 



time constant for out-scattering found from the quantum Boltzmann equation ( 31 1. Moreover, in 
the low density limit, the phase build-up will always be negligible compared to the phase decay 
on average, because the build-up rate (|56| is proportional to the square of the occupation number, 



while the decay term ( 66 1 is linear in the occupation number at low density. 

Turning to the p (4 'terms, which we found can also contribute to the evolution of {Nk), it 
should not be hard to see that the same procedure applied to these phase factors will give the 
same behavior. In the above calculation for the dephasing, we could have just as easily inserted 
p (4) , with the same result, but with eight scattering terms instead of four. In the same way, a 



first-order calculation like the one which gave us the result ( 56 1 for the accumulation of phase, 
presumably will also give us a term linear in t, but this phase accumulation will quickly be 
removed by the dephasing. 

We can therefore say in general terms that the natural evolution of the system makes the 
many-body states tend rapidly toward becoming diagonal states, which was the assumption used 



implicitly in the iterative solution of the quantum Boltzmann equation. In equilibrium, (56 1 
implies that there is no build up of phase factors, and the states resolve to diagonal states, that 
is, states in which only (Nk) is relevant. If the system is far from equilibrium, there will be a 
tiny buildup of phase in the superposition of states which is constantly being destroyed by the 
out-scattering processes. 

Again, note that this entire calculation has been done for a closed system with energy con- 
servation. No notions of measurement or collapse have been invoked. Rather, the deterministic 
evolution of the quantum many-body wave function in a large system gives rapid dephasing 



which leads to information loss. The strongest assumption we have made is that of Section 2.2 
that there is no correlation of the occupation numbers in different Fock states in an overall super- 
position. A similar assumption is made in Section [33] that there is no correlation of phase and 
number in different states. 

Is the system reversible, then? Another way to put this question is whether we could choose 
some special initial state that would cause the system to evolve away from equilibrium and toward 
our initial state. The answer is yes, but our calculation gives an important qualification. To make 
the system evolve backwards, i.e., away from equilibrium, we would have to choose a special 
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set of phase coherences in a superposition of states of the form Q that corresponded to the 
state of the system after some time evolution, set the system to have exactly these phase factors, 
and reverse the time evolution. But these phase factors tend exponentially quickly toward zero. 
Picking the right phase factors therefore is analogous to picking the proper initial state in a 
classical chaotic system to make it evolve toward a predetermined goal. As with predicting the 
motion of a classical chaotic system, time-reversing an interacting quantum many-body system 
involves exponential sensitivity to initial conditions (known as "stiffness") 1 19|. 

In all of this discussion, we have used the case of two-body, elastic interactions. It is not hard 
to see that the main results here do not depend on the exact form of the interactions. Any closed 
many-body system, even one in which particle number is not conserved, will evolve toward 
equilibrium and will have rapid dephasing of the information needed to reverse its evolution. 

The above perturbative approach, in which we limit the calculation to the second order of 
perturbation theory, also assumes weakly interacting particles. If we did not assume weak inter- 
actions, then we could not treat the single-particle states as nearly eigenstates, as we have, and 
talk of particles scattering from one state to another. In other words, the whole particle picture 
would break down. In this case, what is typically done is to renormalize the system to define 
new eigenstates, which then become the new particles in the system (known as "quasiparticles"), 
which have weak interactions between them. In this approach, a strongly interacting system can 
be mapped to a weakly interacting system with a new type of particle. (For example, plasmons 
are the new type of weakly interacting particle in a gas of strongly interacting electrons.) If this 
is not possible in some systems, then we cannot speak of equilibration of particles at all, because 
particles are no longer the proper description of the system. 

This calculation tells us something more than the classical Boltzmann scattering integral with 
the Stosszahlansatz. The quantum approach here says, first, that if we know that the system is in 
a purely diagonal state, then it will evolve deterministically toward equilibrium, and second, that 
the system naturally evolves toward diagonal states on time scales comparable to the scattering 
time. After starting in a nonequilibrium diagonal state, at a later time a system will not be 
in exactly a diagonal state; otherwise we could not time-reverse its state, which would violate 
the underlying time-reversibility. If we do not have exact knowledge of the initial state of the 
system, we cannot rule out the possibility that the initial state of the system was a state with all 
the phases exactly such that it will evolve away from equilibrium. As with the Stosszahlansatz 
in the classical H-theorem, we assume it is extremely unlikely to have been in such a state. But 
since we know such states are exponentially suppressed in time, our choice of a diagonal state as 
the initial state in the absence of complete knowledge of the phases is not arbitrary, but naturally 
favored by the system itself. 

6. Conclusions 

The picture that arises from quantum field theory is that the H-theorem and the Second Law 
arise, not from statistical arguments, but simply from the coupling of local energy to the whole 
vast universe. If a system could truly be isolated perfectly, as a finite, energy-conserving system, 
then it would indeed undergo periodic motion forever. As soon as one allows any coupling 
to the rest of the universe, however, then there is damping which leads things to run down. 
This result has been deduced using exclusively the wave picture of quantum mechanics, without 
reference to wave-particle duality or the stochastic nature of identifying particle events with the 
wave function. No random numbers were used in the numerical calculations. The H-theorem 
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fundamentally comes from the nature of wave behavior in a system with an infinite number of 
degrees of freedom. 

As discussed in the introduction, work by Anatoli Polkovnikov [6 1 has shown that a "diagonal 
entropy," or d-entropy, can be defined by summing only over diagonal states, and this entropy is 
extremely useful in calculations of thermodynamics properties in quantum systems. The present 
work shows that this choice of d-entropy is natural because the off-diagonal terms vanish quickly 
due to dephasing. The present work also bears some similarity to the work of Zurek 11201 and 
others which shows that quantum states can dephase to diagonal states due to interaction with an 
environment; in this work the "environment" is part of a closed, energy-conserving system, and 
all the dephasing is internal. 

Nuclear magnetic resonance spectroscopy provides an example of how this works. Suppose 
that a pick-up coil detects the global ^-component of the magnetization (corresponding to a 
particular superposition of up and down spin) created by an initial n/2 tipping pulse. Initially 
after the tipping pulse, the spin up and spin down parts of the wave function have the same 
spatial function. However the spin up and spin down components of the particle take different 
trajectories and so eventually spatially separate from each other, entangling the spin with orbital 
motion. The spin-up part of one packet can collide with the spin-down part of another packet, and 
these collisions spread the entanglement throughout the system just like the collisions we have 
considered in this paper. The time evolution is completely unitary but produces multi-particle 
entangled states which look locally (i.e., in terms of individual spins) to be incoherent mixed 
states. That the evolution is unitary and ultimately coherent can be seen in sophisticated echo 
techniques which can be used to narrow the resonance lines by orders of magnitude and recover 
the signal ||2T|| . 

Another example can be found in ultra-cold atomic gases which are quite isolated and essen- 
tially uncoupled from any heat baths. Here the natural coherence is so high that it can actually 
be experimentally difficult to produce incoherent mixtures of different spin states. In the case 
of cold fermions with short-range interactions, the Pauli principle renders these interactions in- 
effective except between opposite-spin particles. Starting with all spin-up particles, application 
of a uniform n/2 pulse on all the spins does not turn on the interactions, even though both spin 
up and spin down are now present, because all the particles are in the same coherent superpo- 
sition spin state. However, in the presence of a spatial gradient in the Zeeman field, the spin 
state of the moving particles eventually becomes entangled in a complicated way with the spatial 
orbital state. Despite the overall unitary time evolution, this leads to an effective decoherence 
of the spins, which turns on the interactions and leads to observable mean-field shifts in radio 
frequency spectroscopy experiments [22 23l l24l . 

The above discussion in terms of off-diagonal terms of the generalized density matrix is 
made a little more complicated by the phenomenon of spontaneous phase coherence, which is 
seen experimentally in boson systems such as superfluids, superconductors, and lasers. In such 
systems, as a system evolves toward equilibrium, a macroscopic fraction of the whole system 
can obtain a finite off-diagonal density matrix term (see Appendix C). This fraction is known as 
the "condensate." This leads to periodic behavior on long time scales, seen for example in the 
fascinating effects of persistent current and Josephson oscillations of superconductors. 

Furthermore, all classical wave states such as water waves on the ocean or radio waves can 
be viewed as condensates with macroscopic phase coherence. This is because all classical waves 
are technically macroscopic coherent states of bosons. (See Ref. [9], Section 4.4, and Ref. 11251 . 
for discussion of coherent states, also known as Glauber states.) The symmetry breaking in 
this case does not come from a thermodynamic phase transition, as occurs for superfluids and 
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superconductors, nor even from amplification of small fluctuations, as in a laser, but from an 
external driving force; thus classical waves are sometimes called "driven condensates" 11261 . The 
external driving force breaks the symmetry and favors one definite phase over all others. If 
the driving force is removed, however, the interaction of many degrees of freedom eventually 
leads to dephasing. Classical waves typically can be sustained because the bosons of which they 
are comprised have very weak interactions; the phonons in sound waves have extremely weak 
interactions at low frequency (see Ref. (9), Sections 5.3 and 5.4), while photons in vacuum are 
essentially non-interacting. This allows coherent states to be sustained for long periods of time. 
Eventually, however, in the absence of an external driving force, they dephase and turn into a 
diagonal Planck equilibrium distribution, which we call heat. 

Condensates and classical waves have, to a high degree of approximation, reversible behavior. 
No small system can ever be perfectly isolated from the rest of the universe, however. Because 
of tunneling, there will always be coupling to the infinite external universe, which will show 
up as damping in smaller systems. Thus all local systems dephase eventually due to radiation 
loss to the infinite space of the universe. This may occur extremely slowly, as in, for example, 
the cosmic background radiation which has not yet perfectly equilibrated over the eons, but the 
Second Law still holds. 

This formulation of the Second Law has important cosmological implications. A prevalent 
model of the universe to account for the data of accelerating expansion [27] is that it is open, 
which means that it is infinite in spatial extent on any equal-time slice [28 1, essentially an infinite 
ensemble of causally disconnected spaces. The possibility of an truly infinite ensemble can 
lead to bizarre implications if the standard, statistical version of the Second Law is used. For 
example, the fact that stars radiate energy rather than sucking energy in is due to the Second 
Law. If the universe is truly comprised of an infinite ensemble of locations, then on purely 
statistical grounds we would expect all finitely improbable events to occur, not just once, but 
an infinite number of times (though not necessarily observable to us). Thus we might expect, 
for example, that somewhere in the universe there are stars that suck in light rather than emit it 
(for any finite length of time). The quantum mechanical formulation presented here says that if 
stars are coupled to an infinite universe, then the Second Law is absolute, not just a likelihood, 
unless the beginning of the universe had very bizarre initial conditions of phase coherence. Any 
finite system coupled to the infinite background of the universe will have the same deterministic 
irreversibility. 
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Appendix A. Connection of the Boltzmann Equation to the H-Theorem and Entropy 

Here we derive the H-theorem from the quantum Boltzmann equation, following the standard 
approach based on the classical Boltzmann integral, e.g. Ref. PP . 

For a quantum system, the total entropy is given by the von Neumann entropy 

S vN = -k B lxip\np\ (A.l) 

where p is the full density matrix of the system. This is not the same as the single-particle density 



matrix defined in (44 1, but instead is defined for the set of all Fock states. Using the general state 
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we have p mn = a* n a m . Since the quantum mechanical evolution of a closed system involves 
only unitary operators, this entropy never changes. 

Polkovnikov |6| has suggested another entropy quantity, which can be called the diagonal 
entropy, 

S d = -k B J](N k }ln(N k }. (A.2) 

k 

This corresponds to the standard thermodynamic entropy for an ensemble of classical particles. 
Like the concept of entanglement entropy 0, it involves a subset of the full quantum mechanical 
wave function information which is relevant to the macroscopic physical properties of the system. 
It also has relation to the concept of quantum coarse graining [5|, because it involves averaging 
over many states. 

The time derivative of the diagonal entropy ( A.2 1 is 

The second term vanishes due to number conservation. The first term is given by the quantum 
Boltzmann equation ( 3 1 1, so that we have 

— - = -k B C J] ln<i9*> [<K )(N k2 )d ± <&»(1 ± 

-(N k )(N k ,)(l ± (N ki m ± <^ fa »] ■ (A.4) 

where C is a factor that contains the matrix element and conservation of energy factor. 

If we pick four particular states k, k' , k\ , and k%, then the total of all terms in the sum involving 
the same four states is 

(ln<# t > + ln<#*> - ln<^,) - ln<# fa >) 

x [0 h )(ft k2 )(l ± (N k ))(l ± - (N k )(N k ,)(l ± <J9fe»(l ± <^ fe »] 

= ln(f^'>) (A.5) 
x [(ft h )($h)a ± <^»d ± <^» - (N k )(N k ,)(l ± <J9fe»(l ± <# fe »] . 
If the density is low, so that the ((1 ± (N k )) factors are negligible, this becomes 

ln ( < *>< >> \0 ki )( # ■ y _ { N k}{ N k , } ] , (A.6) 

which is always negative or zero for values of the occupation numbers greater than zero. Since 
the whole sum (A.4i consists of terms like this, the total sum is less than or equal to zero, and 
therefore dSd/dt > 0. This is the standard form of the H-theorem. 

If the system obeys fermion statistics, the full sum ( A.5 1 is also always negative or zero. But 
if the system obeys boson statistics which allow occupation numbers greater than unity, it can be 
made positive; for example, if (A^) = 10, (N k >) = 0.1, (A*,) = 0.5, and {N k2 } = 0.5. Therefore, 
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for some distribution functions dSd/dt < 0. To avoid this problem, the proper semiclassical en- 
tropy quantity analogous to \A.2\ should be used. From standard quantum statistical mechanics 
(32), this is 

Ssc = -kB^imMm + ii + iminii + iNk))), (a.7) 

where the upper sign is for bosons and the lower sign is for fermions. Assuming number conser- 
vation, the time derivative of this is 




(A.8) 



Using the quantum Boltzmann equation for d{Nk)ldt as above gives dS sc /dt > for all distribu- 
tions. 



Appendix B. Calculation of the Off-Diagonal First-Order Terms 

the first-order evolution of the p (2) operator gives us terms of the 



As discussed in Section 
following form: 



5.2 



{\p (2) (t), Vi(t')]) = UinV^alala^e-^e^al a\ a k a k e- iH °<'\n) 



i k u k ,u k „u k „,e e u k u k u k2 u h 
2 k 4 CL k, a k 2 a k l e e Q k a k>' 



-U(n'\e iH <"'al a\ a k a k e~ m ° f ^alal.a^a^e'^n). (B.l) 



We must treat a number of different possibilities: 1) one of the creation or destruction operators 
in Vj does not commute with one of the operators in p (2) ; 2) two of the operators in V, do not 
commute with p (2) ; 3) three of the operators in V, do not commute with p (T> ; or 4) four of the 
operators in V,- do not commute with p <2 \ We take these separately; each has several possibilities. 

One non-commuting operator. If there is one non-commuting operator in V,, we will have 
terms like the following: 

{a\a k - ^t fl ^^^^^-^W e -« + ^-^/» 

= -a\,alala k ,,a kl ,,a h e i{ ^ (B.2) 

This is proportional to a third-order off-diagonal term fP> with six creation and destruction opera- 
tors. We ignore these higher-order off-diagonal terms for the moment; in the spirit of perturbation 
theory, we can expand in a series of orders of off-diagonal terms, with the most significant being 
p^, assuming that as the degree of correlation increases, the lower the magnitude of the term. 

Two non-commuting operators. There are three cases: a) two operators make up a number 
operator N k ; b) two operators act on different fc-states; or c) two operators act on the same state 
and are the same. 

a) In this case, we will have terms of like the following: 

{a[N k -N k a\)a[,a kl ,a k ,,,a[^ 
= -^ v a v ^^^ E ^- E ^e^- E ^ l ' ,l . (B.3) 

If K2 + k?,, we again have a higher-order p <3) term. If ki - £3, then this becomes 

- a[a[,a k „a k ,,,N k J (E ^ Ek '- Ek "- Ek, " )t!r \ (B.4) 
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But in this case there is always also another term 

{a k „A»' - fitra^yfykarfik^**-^- 1 ^* 1 * (B.5) 

which cancels the first term. This is true for both bosons and fermions. 
b) In this case, we will have three types of terms. The first is: 

(a i aa t a - a a f a aha a a f J(E k +E k ,-E t ,,-E k „,)t/ti -i(E k +E k ,-E h -E M )f/n 
\a,a k a.,a k , a k ,a k ,a k a,)a k „a k ,„a, a. e e 



(B.6) 



The time-dependent exponential factors in lowest order become just t by the same type of calcu- 
lation as in ( |55| ). So we have 

- (+N k ± fa + \)a kll a k „ ia \a\^ (B.7) 

where the ± signs give the + sign for bosons or - sign for fermions. 
The second type is: 



{a k „a\„a k ,„a\,„ - a\,,a k ,,a\,„a k ,,?,a\a\,a k a k t = {±N V , ± N k „ + \)a\a\,a k a k t, 



(B.8) 



and the third type is: 



{a\a k a K „a\„ - a k a\a\,,,a k ,^ k ,a k ,,a^ k a ki t = (N k - N k »<)a} k ,a k „a\a k t. (B.9) 

All of these depend on various p (2) factors. In general, they do not cancel, but one can expect 
that they average to zero when the sum over the free momenta is done. 

c) These terms vanish for fermions and give higher order p (3) terms for bosons. 

Three non-commuting operators. There are two possibilities here: a) all three act on dif- 
ferent fc-states, b) two act on one state, and the other acts on a different state. 

a) In this case, we will be left with factors that depend on the occupation numbers N k times a 
two-operator off-diagonal element of the form p k kl = aia k . If the system is in a diagonal state, 
this two-operator off-diagonal element cannot become nonzero through the two-body interaction 
we consider here. However, if there are nonzero p (2) terms, then the first-order evolution of 
can be nonzero. 

b) This will lead to terms which depend on various p (2) factors. 

Four non-commuting operators. This is the same as the case when all the states in \n) are 



restored, treated in Equation ( 56 1. 



Appendix C. Onset of Phase Coherence in a Bose-Einstein Condensate 

The rate of onset of phase coherence in a homogeneous condensate can be found be a simple 



extension of the second-order rate equations for phase already presented in Section 5.3 Instead 
of looking at p^ 2 ], ,,, ,,„ = a k a k ,a k „a k ,„, we instead just insert a k in the same equations. In the case 
of a condensate we have just one term, 

j f {a k ) = (a k )^(U D + U E f l - J] 6(E k + E kl _ - E h - E h ) (C.l) 

kite 

x [(N h )(N k4 )(l + <# fa » - <# fa >(l + (N h W + <#*,»] ■ 
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When there is net influx into the condensate state, any small coherence will be amplified, with 
expontential growth of the amplitude of the phase. The growth of the coherent amplitude will 
cease in equilibrium, whenA^ + l - N k . Then the condition N^N^il+N^) = A^^, ( 1 + A^ 3 )( 1 +A^ 4 ) 
is equivalent to N k ,Nk 4 (l + Nk 2 )(l + Nk) = N k N kl (l + N k} )(l + iVfc,), which is equal to zero in 
equilibrium. 

The calculations for p (2) can also be used to find the evolution of the phase of a condensate. 



The zero-order term ( 5 1 1 corresponds to 



d{a k ) = -i — t(a k ) = -iu k t{a k ). (C.2) 
n 

In other words, the phase of the condensate rotates in time. The first-order term ( |53] l corresponds 
to 



d(ipi\a k \if/i) 

m 



X ({n'\a k e iH ^^/ k2 a h a ki e- iH ^ n \n) - (n'\e iH °' ' l n a\a\a k ,a u e- m ^ '' n a k \n)) . 



The terms with k\ -k^-k and kq, = k$, plus the equal terms with ki—k->,-k and k\ = £4, give 
us 



d{^i\a k \\pi) 



1 r' 

= I dt'2^a* n ,a„2^ Ue~ lEk,,h (n'\{a k a\ - ala k )N h a k \n). 

" n,ri k 2 

- -S-UNMaMi), (C.4) 
n 



where in the last step we took the limit e~' Ek ' /n ^ 1, for short times. This term gives us the first- 
order, mean-field renormalization of the condensate energy, which gives a shift of the frequency 
a>k of rotation in ( |C.2| >. 

The time evolution which leads to onset of phase coherence in condensates has been an 
active topic of study over the past two decades ||29l . The results here apply to a homogeneous 
gas in the thermodynamic limit, showing that spontaneous symmetry breaking occurs in the 
thermodynamic limit on short time scales, not only in the finite systems studied experimentally. 

Several decades ago, Uhlenbeck lf30ll pondered the problem of using the assumption of coarse 
graining to justifying the Second Law of thermodynamics in connection with the macroscopic 
wavelike behavior of condensates. His intuition was basically correct: the assumptions which 
lead to irreversible behavior in normal systems break down in the case of condensates. The role 
of the coarse-graining assumption in the thinking of Uhlenbeck's day is here seen to be played 
by the fast-dephasing assumption. 

We note that the above type of calculation applies to any Bose system, including ones without 
number conservation, such as phonons which a Planck distribution in equilibrium. Whenever 
occupation numbers exceed unity, there can be amplification of coherent phase fluctuations. It 
can be argued that this is the fundamental reason why macroscopic systems such as water waves 
are always in definite-amplitude states instead of Fock states. A macroscopic Fock state, which 
is a superposition of all different phases, i.e. a water wave which simultaneously has a crest and a 
trough, is a physically possible state. The reason such is never seen is related to the spontaneous 
increase of phase coherence in boson systems. 
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